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ABSTRACT 

The narrow peak recently found in various pionic double charge exchange (DCX) cross 
sections can be explained by the assumption of a universal resonance at 2065 MeV, called 
d! . We calculate the mass of a six-quark system with J p = 0~, T = quantum numbers 
employing a cluster model and a shell model basis to diagonalize the nonrelativistic quark 
model Hamiltonian. 
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1. Introduction 

Recently, a very narrow peak (T mediurn rs 5 MeV) has been observed in the ) double 

charge exchange cross section at an incident pion energy of = 50 MeV and forward angles 
e = 5° 0. Dedicated DCX experiments on a number of nuclei ranging from 12 C to 56 Fe have 
shown that the position and width of this peak is largely independent of the nuclear target. 
Therefore, a fundamental two-nucleon or six-quark process seems to be involved. In fact, due 
to charge conservation the DCX reaction involves at least two nucleons within the nucleus 
and is quite sensitive to short-range NN correlations 0. Conventional DCX calculations 
have so far been unable to explain these experimental results I. On the other hand, the 
assumption of a single narrow resonance, called d! , with spin, parity J = 0~ and isospin 
T = and a resonance energy of M# = 2065 MeV works extremely well in describing all 
available DCX data 0. The mass of the cf-resonance coincides with = 2M^ + m n + T n Q. 
The small decay width of the d' is naturally explained by its small mass which is below any 
baryon-baryon (NN* , N*N* , ...) threshold and by its quantum numbers J p =0~ and T = 
which prevent a decay into the nucleon-nucleon channel. 

QCD does not exclude a rearrangement of the six quarks into more exotic configura- 
tions such as a diquark (g 2 ) and a tetraquark (g 4 ) cluster. In fact, the bag-string model 
predicts that a diquark cluster with spin (isospin) S = (T = 0) and a tetraquark cluster 
with spin (isospin) S — 1 (T — 0) moving with relative angular momentum L — 1 is en- 
ergetically the most favorable configuration for a system with d! quantum numbers. Ref.0 
predicts a d! mass of ~ 2100 MeV. However, in these calculations the quark exchange 
interactions between the clusters have been neglected. Therefore, we perform a calculation 

*Tt has been suggested that the d'-dibaryon is a tyNN resonance with isospin T = 2 0. However, the T — 2 
assignment seems to be in conflict with experimental information from other sources 0. 



in the nonrelativistic quark model (NRQM) in which the complications due to the Pauli 
principle are properly taken into account. We use two different sets of basis states (i) a clus- 
ter model basis using the Resonating Group Method (RGM), (ii) a quark shell-model basis. 
Furthermore, unlike the bag-string model, the NRQM calculation includes the nonperturba- 
tive ii- and a-meson exchange interactions between quarks as required by the spontaneously 
broken chiral symmetry of QCD. 



2. The nonrelativistic quark potential model 

In the NRQM a system of n-quarks with equal masses m q = 313 MeV is described by the 
Hamiltonian 



H = E(m q + ^)- + E V C ° nf (^ r,) + £ V»(r 4 , r,), (1) 

where rj, p« are the spatial and momentum coordinates of the i-th quark, respectively and 
P is the total momentum of the n-quark system. The exact removal of the kinetic energy of 
the center of mass motion by the third term is crucial. 

The long-range confining force is generated by a two-body harmonic oscillator con- 
finement potential 

V conf (r h r j ) = -a c X i -X j (r i -r j ) 2 , (2) 

where Aj is the Gell-Mann matrix of SU(3) co i or of the i-th quark. The residual interaction 
models the relevant properties of QCD, namely asymptotic freedom at short distances and 
spontaneous chiral symmetry breaking which has important consequences for the short- and 
intermediate-range interaction between quarks. The one-gluon exchange potential 

V OGEP ^ rj) = f A, • A,{i - ^ (l + \a, • a,) *(r)} (3) 

provides an effective quark-quark interaction that has the spin-color structure of QCD at 
short distances. Here, r = r 4 — and o"j is the usual Pauli spin matrix. 

It is well known that the spontaneous breaking of chiral symmetry by the physical 
vacuum is responsible for the constituent quark mass generation i, as well as for the ap- 
pearence of almost massless pseudoscalar Goldstone bosons and massive scalar mesons that 
couple to the constituent quarks. In the NRQM, this is modeled by regularized one-pion 
and one-sigma exchange potentials between constituent quarks 0: 
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|f = |f, ml = {2m q f + ml K = A a = A. (5) 
The nq coupling constant is simply related to the tcN coupling constant H' 0. We use 



f 2 N /47i = 0.07491111. The nq cut-off mass A describes the extended pion-quark vertex. Here, 



A = 4.2 fm 1 which results in a soft ttN form factor 113. 



We take a harmonic oscillator wave function for the i-th quark inside the nucleon and 
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where b is the harmonic oscillator constant. 



Table 1: Quark model parameters. Set I: with regularized 7r- and cx-meson exchange poten- 
tials. The a parameters are fixed by eq.([|). Set II: without ir- and a-exchange potentials 
111. Set III: as set I but with reduced confinement strength (see sect. 4). 



Set 


b [fm] 


a s 


a c [MeVfm- 2 ] 


m a [MeV] 


£/(4f) 


A [fm- 1 ] 


I 


0.595 


0.958 


13.66 


626 


0.554 


4.2 


II 


0.603 


1.540 


24.94 








III 


0.595 


0.958 


5.00 


626 


0.554 


4.2 



We determine the parameters a c , a s , and b by requiring that the AT(939) and A(1232) masses 
are reproduced and that the nucleon mass is stable with respect to variations in b. 



M N {b) = 3m q = 939MeV, M A - M N = 293MeV, 
3. Six-Quark Models for the d'-Dibaryon 



dM N (b) 
db 



0. 



(7) 



In the bag-string model, dibaryons are described as rotating strings with colored quark 
clusters at the ends i. This model leads to a linear Regge trajectory of excited states. It 
predicts that the lowest L = 1 excited state is obtained for a diquark-tetraquark configuration 
at around 2100 MeV I. 




Figure 1: The d'-dibaryon in the q -q quark cluster model. 



The drawback of this approach is that it does not respect the Pauli principle. Only the 
quarks within the individual clusters are antisymmetrized but not the quarks belonging to 



different clusters. This is a good approximation for high angular momentum states because 
in this case the system is fairly elongated and the probability of cluster overlap is small. On 
the other hand, for a low lying L = 1 excitation, such as the d! one expects a considerable 
amount of quark exchange between the two clusters. The confining forces between the 

Table 2: The mass (M<f) and size (6 6 ) of the d! in the quark cluster model. The masses of 
the diquark and tetraquark are also given. 





M 2q 


M 4q 


M d , 


h 


Set I 


645 


1455 


2440 


0.76 


Set II 


637 


1501 


2634 


0.70 


Set III 


621 


1309 


2111 


0.95 



colored quarks prevent large separations of the clusters and the typical size of such a system 
is expected to be about 1 fm. From our experience with the iViV system we know that the 
Pauli principle plays an important role at such distances 0. 

3.1 The Quark Cluster Model of the d'-Dibaryon 

In this model, the d! is described as a nonrelativistic six-quark system in which the quarks 
interact via the two-body potentials of eq.(^). Tensor and spin-orbit interactions have been 
omitted since it has previously been shown that they give a negligible contribution to the dl 
mass El. The six-quark wave function is expanded into the cluster basis 

| ^=0,T=0 > = A 

^= o ' T - o (p D )xgj 5 ^°® Xi=1 (R) 

where ^ T ~ 1,TT ~°(p T , \ T ,rf T ) and $b d ~°' Td ~°(Pd) are the internal wave functions of the 
tetraquark (T) and diquark (D) clusters, respectively and Xl=i(R-) is the wave function of 
the relative motion of the two clusters. We use the same harmonic oscillator parameter 
for the internal and relative motion wave functions. The Young diagrams in eq.(G|) show 
that two color triplet clusters [211]c and [ll]c are coupled to a [222] c color-singlet six- 
quark state. Furthermore, they show that the tetraquark and d! wave function are not fully 
antisymmetric but have mixed symmetry in color space. This nonfactorizability of the color 
space considerably complicates the calculation. 

The advantage of the cluster model is that it provides a continuous transition from the 
q 6 six-quark state to the q 4 — q 2 clusterized state by smoothly going through all intermediate 
configurations. There is no rigid and artificial boundary between these extreme configura- 
tions; everything is contained in one and the same RGM wave function. This important 
property is a direct consequence of the Pauli principle on the quark level, which is ensured 
by the antisymmetrizer A 
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Figure 2: The direct, one-quark and two-quark exchange diagrams that have to be evaluated 
for each two-body potential. The horizontal bars indicate the confinement, the one-gluon, 
one-pion, or one-sigma exchange interactions in eq.( [I]). 



where P® STC is the permutation operator of the i-th and j-th quark in orbital (O) spin- 
isospin (ST) and color space (C), respectively. The direct, as well as the one- and two-quark 
exchange contributions for the two-body potential of eq.([l]) are depicted in fig. 2. The solution 
for the unknown relative wave function Xl(R) an d the unknown eigenenergy is obtained from 
the variational principle 

-i^ d ,\H-E\^ d ,y 
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(10) 



where the variation is with respect to the relative wave function x^(R). The results for the 
energy (mass) of the d! as well as for the harmonic oscillator parameter 6 6 which minimizes 
the d' mass are shown in table 2 for the parameter sets of table 1. 



3.2 Shell-Model Calculation for a J p =0 , T=0 six-quark system 

Next, we calculate the mass of the ci'-dibaryon in the translationally invariant shell- model 
(TISM) II. Due to the negative parity of the d', only an odd number of oscillator 
quanta N = 1,3,5, ... is allowed. There is only one N = 1 state which is compatible with 
jp = 0",T = 

| t> dlg s >=\N = 1, [51] , (A/i) = (10), L = 1, S = 1, T = 0, [321] ST > . (11) 

For an unambigious classification of TISM states one has to specify the number of internal 
excitation quanta N, the Elliot symbol (A/i), the Young pattern [f]o of the spatial permu- 
tational S^-symmetry, further the total orbital angular momentum L, total spin S and total 
isospin T of the system. The specification of the intermediate SU(4:)st symmetry is nec- 
essary because in general, the same symmetry in STC space can be obtained from several 
states with different intermediate ST symmetries. The mass of the d' is then given in first 



order perturbation theory by the expectation value of the Hamiltonian between the lowest 
harmonic oscillator state of eq.flTll) 

M d ,{h) =< % s | H | * dlga > . (12) 

In order to estimate the effect of configuration mixing with excited shell model states we 
include in addition ten N = 3 states with orbital [42] o symmetry &3. In this case also the 
[51] st, [411]st, [33] st, [321]st, and [2211]st S§ permutational symmetries are allowed. 

With fixed parameters of the quark-quark interaction determined from eq.(|7|) we 
minimize the d' mass with respect to the harmonic oscillator parameter 6 6 in the six-quark 
wave function. Note, that the harmonic oscillator parameter of the single baryon (b) and the 
d! (bo) wave function are different. The value of b^ which minimizes the d' mass is a measure 
of the size of the system and is also given in table 3. 

Table 3: The mass (M d >) and size (bo) of the d' in the six-quark shell model without and 
with configuration mixing. 





N = 1 


N= 1 & N = 3 




M ( f =1) [MeV] 


h [fm] 


M d . [MeV] 


h [fm] 


Set I 


2484 


0.78 


2413 


0.78 


Set II 


2636 


0.72 


2553 


0.73 


Set III 


2112 


0.95 


2063 


0.96 



4. Discussion and Summary 

As is evident from tables 2 and 3 the d! mass of the cluster model is lower than the single 
N = 1 shell-model mass of eq.(^) but higher than the shell- model result with configuration 
mixing. Note that the d' mass calculated with set II (without pion and sigma-exchange 
between quarks) is some 150-200 MeV higher than the result with chiral interactions (set 
I). In any case, the calculated mass is about 350 MeV higher than the value required by 
experiment. However, the confinement strength a c in the three-quark and six-quark system 
need not be the same. If we assume (set III) that a c in the six-quark system is weaker than in 
the nucleon one obtains considerably smaller values for M d >. This assumption is supported 
by the harmonic oscillator relation for a c 

« C oc — -— (13) 

which is inversely proportional to the number of quarks iV in the system. A weaker con- 
finement strength is also expected due to the larger hadronic size of the d! ib§) as compared 
to the hadronic size of the nucleon (b). Set III differs from Set I of table 1 only in the 
strength of the parameter a c for which we take the value a c = 5.0 MeV/fm 2 in the six-quark 
calculation. Finally, both calculations give similar results for the d! mass and for its size. 
Let us briefly discuss the reasons for this. The outer product of the orbital [4]o (tetraquark) 
and [2]o (diquark) symmetries gives the following six-quark symmetries 



[4]o ® [2]o = [42] © [51]o © [6]o- 



(14) 



With the exception of the [6]o symmetry which is incompatible with d 1 quantum numbers 
these are also included in the enlarged N = 3 shell-model basis H. Analogously, the outer 
product of the two clusters in spin-isospin space leads to 

[31] 5T ® [2)st = [5l) ST © [42] 5 t © [33] 5 t © [411] OT © [321]^. (15) 

Comparison with eq.(10) in refS shows that the g 4 — q 2 cluster model wave function com- 
prises the same ^symmetries in orbital and spin-isospin space (with the exception of the 
[2211]5r symmetry) as our enlarged shell model basis. Thus the trial function space spanned 
by both sets of basis functions is not very different. 

In summary, we have calculated the mass ofaJ p = 0~T = six-quark system in the 
NRQM using two different assumptions for the spatial distribution of the six quarks. The 
parameters have been determined from the constraints of eq.(^). As in our previous works 
B 11 our results are typically 300-400 MeV above the required resonance energy. However, 
for a weaker confinement strength a c in the six-quark system as suggested by eq.(|T3D we 
find a mass for the d' that is considerably smaller. The assumption of a weaker confinement 
strength in the six-quark system does not affect previous results of the model in the B = 2 
sector such as NN scattering phase shifts or deuteron electromagnetic form factors which 
are completely insensitive to the model and strength of confinement 0. 
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